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J  -  j  Abstract  M — LA - 1. 

We  review.some  recent  mathematical  results  concerning  integrodifierential  equations 
that,model  the  motion  of  one-dimensional  nonlinear  viscoelastic  materials.  In  particular, 

v  .  •'  J 

we  discuss  global  (in  time)  existence  and  long-time  behavior  of  classical  solutions,  as  well 

as  the  formation  of  singularities  in  finite  time  from  smooth  initial  data.  Although  the 
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mathematical  theory  is  comparatively  incomplete,  jyenmak^some  remarks^concerning  the 

existence  of  weak  solutions  (i.e.,  solutions  with  shocks).  Some  relevant  results  from  linear 

Is  ,  .  r  *  _  .'/ ..  ^  > 

wave  propagation  will  also  be  discussed.  rf  r 
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INITIAL  VALUE  PROBLEMS  IN  VISCOELASTICITY 
W.  J.  Hrusa’,  J.  A.  Nohel",  and  M.  Renardy’" 


1.  Introduction 

Over  the  last  decade  or  so,  a  significant  amount  of  research  has  been  devoted  to  the 
study  of  integrodifferential  equations  which  model  the  motion  of  nonlinear  viscoelastic 
materials.  The  purpose  of  this  paper  is  to  review  some  of  the  results  that  have  been 
obtained  and  some  problems  that  remain  open.  We  focus  on  equations  that  describe 
viscoelastic  solids;  a  similar  paper  concerning  viscoelastic  fluids  will  appear  elsewhere. 

Our  presentation  will  be  informal  in  the  sense  that  technical  hypotheses  of  smoothness, 
etc.  will  not  always  be  stated  explicitly,  and  we  shall  not  attempt  to  state  results  in  the 
greatest  possible  generality.  We  refer  to  the  original  sources  and  to  our  recent  monograph 
(Renardy,  Hrusa,  and  Nohel  (1987),  henceforth  referred  to  as  RHN)  for  precise  statements 
of  technical  assumptions  and  for  more  general  results.  In  order  to  avoid  geometrical  com¬ 
plications,  and  because  the  mathematical  theory  of  such  equations  is  far  more  complete, 
we  restrict  our  attention  to  problems  involving  a  single  spatial  coordinate.  Moreover,  we 
always  treat  the  full  dynamical  equations  with  inertia;  quasistatic  approximations  will  not 
be  discussed. 

Consider  the  longitudinal  motion  of  a  homogeneous  one-dimensional  body,  e.g.  a 
bar  of  uniform  cross-section.  [This  problem  is  somewhat  artificial  since  a  longitudinally 
deforming  bar  may  experience  a  variation  in  its  cross-section.  There  are,  however,  more 
realistic  classes  of  motions  that  can  be  described  with  a  single  spatial  coordinate  (cf.,  e.g., 
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the  appendix  of  Coleman,  Gurtin,  and  Herrera  (i9G5)  and  Section  1.4  of  RHN).]  Wc  assume 
that  the  body  occupies  an  interval  B  C  III  in  an  unstressed  state.  The  interval  B  will  be 
called  the  reference  configuration;  a  typical  point  (or  particle)  in  B  will  be  denoted  by  x. 
In  order  to  describe  the  motion,  we  follow  the  evolution  of  points  in  B.  Therefore,  the 
quantities  of  interest  will  be  regarded  as  functions  of  the  material  point  x  and  the  time  t. 

We  assume  that  the  body  has  unit  reference  density,  and  we  denote  by  u[x,t )  the 
displacement  at  time  t  of  the  particle  with  reference  position  x  (i.e.,  x  +  u(x,t)  is  the 
position  at  time  t  of  this  particle).  The  strain  e  is  given  by 

e(z,f)  =  ux(x,t),  (1.1) 

and  the  equation  of  balance  of  linear  momentum  has  the  form 

uu(x,t)  =  ax{x,t)  +  f(x,t),  x£B,  t>  0,  (1.2) 

where  a  is  the  stress  and  /  is  a  (prescribed)  body  force.  Subscripts  x  and  t  indicate  partial 
derivatives;  we  use  a  prime  to  denote  the  derivative  of  a  function  of  one  variable.  The 
type  of  material  composing  the  body  is  characterized  by  a  constitutive  assumption  which 
relates  the  stress  to  the  motion. 

If  the  body  is  elastic  then  the  stress  depends  on  the  strain  through  a  constitutive 
equation  of  the  form 

o{x,t)  =  p(e(x,t))>  (1-3) 

where  is  a  given  smooth  function  with  £>(0)  =  0  and  <p'(0)  >  0.  The  condition  <p(0)  =  0 
reflects  the  assumption  that  reference  configuration  be  stress-free,  while  the  condition 
V?'(0)  >  0  means  that  stress  increases  with  strain,  at  least  near  t  —  0.  For  simplicity,  let 
us  assume  temporarily  that  B  -  IR  and  /  =  0.  The  resulting  motion  is  then  governed  by 
the  quasilinear  wave  equation 


We  seek  a  function  u  that  satisfies  (1.4)  together  with  the  initial  conditions 


u(i,0)  =  uq(i),  u*(x,0)  =  uj(x),  x  €  Ht. 


In  the  linear  case,  a  =  Ee  (with  E  a  positive  constant),  equation  (1.4)  reduces  to 


utt  =  Euxx,  z  6  IR,  t  >  0.  (1.6) 

A  well-known  feature  of  (1.6)  is  that  singularities  in  the  data  (i.e.  discontinuities  in  the 
data  or  in  derivatives  of  the  data)  persist  as  singularities  in  the  solution  for  all  t  >  0 
and  they  propagate  with  constant  speed  \/E.  The  amplitude  of  such  a  singularity  neither 
grows  nor  decays. 

For  the  linear  equation  (1.6),  singularities  will  be  present  in  the  solution  if  and  only 
if  they  are  present  in  the  data.  In  the  nonlinear  case  the  solution  generally  develops 
singularities  (in  finite  time)  -  no  matter  how  smooth  the  data  are  -  because  the  wave 
speed  v^7  is  not  constant.  This  issue  will  be  discussed  in  more  detail  in  Section  4. 

Experience  indicates  that  certain  materials  have  memory,  i.e.  the  stress  depends  not 
only  on  the  strain  at  the  present  time,  but  also  on  the  entire  temporal  history  of  the  strain. 
Typically  the  memory  fades  with  time,  i.e.  disturbances  which  occurred  in  the  distant 
past  have  less  influence  on  the  present  stress  than  those  which  occurred  in  the  more  recent 
past.  We  shall  refer  to  such  materials  as  viscoelastic.  The  linearized  constitutive  relation 
for  small  deformations  was  given  by  Boltzmann  (1876): 

a{x,t)  =  0c(x,t)  +  /  m(f  —  r)(e(x,t)  —  e(x,r))dr.  (1.7) 

J  — OO 

Here  (3  is  a  non-negative  constant,  and  m  is  a  positive  monotone-decreasing  function.  To 
ensure  that  (1.7)  is  meaningful,  we  assume  that  m  is  integrable  at  infinity  and  that  sm(s) 


is  integrable  at  zero,  i.e. 


J  m(s)ds  <  oc,  J  sm(s)d 


s  <  oo. 
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If  m  is  integrate  on  (0,  oo),  i.e. 


fOO 

I  m(s)ds  <  oo, 

Jo 

we  can  rewrite  (1.7)  in  the  alternative  form 

<r(z,t)  =  Ee(x,t)  -  f  m(t  -  T)e(x,r)dT, 
J  —  OO 


(1.9) 


(1.10) 


where 


E  =  3  +  /  m(s)d.s.  (1.11) 

Jo 

(Since  m  is  assumed  to  be  positive,  decreasing,  and  integrable  at  infinity,  (1.9)  holds  if 
and  only  if  m  is  integrable  at  zero.)  The  constant  E  measures  the  instantaneous  response 
of  stress  to  strain  and  is  called  the  instantaneous  stress  modulus.  On  the  other  hand,  3 
determines  the  stress  due  to  a  constant  strain  history  and  is  called  the  equilibrium  stress 
modulus;  for  viscoelastic  solids  it  is  natural  to  assume  3  >  0.  If  we  set  c  =  0  for  t  <  0  and 
e  =  «o  for  t  >  0,  then  the  stress  for  t  >  0  is  given  by 


fOO 

o  =  [3  +  J  m(s)ds)e0  =  G(t)e0. 


(1.12) 


The  function  G  defined  through  (1.12)  is  called  the  stress  relaxation  modulus.  An  experi¬ 
ment  in  which  the  strain  is  suddenly  changed  from  0  and  held  at  a  constant  value  is  called  a 
stress  relaxation  test;  by  measuring  the  stress  one  obtains  G(t).  Obviously,  G(0)  =  E  and 
G(oo)  =  3 ■  The  function  m(t)  =  — G'(t )  is  called  the  memory  function.  The  assumption 
that  m  is  positive  means  that  G  is  monotone-decreasing,  i.e.  that  the  stress  following  a 
sudden  deformation  relaxes  with  increasing  time.  It  follows  from  (1.7)  and  the  positivity 
of  m  that  the  stress  needed  to  sustain  a  strain  c  is  diminished  by  previous  strains  of  the 
same  sign  and  increased  by  previous  strains  of  the  opposite  sign.  The  assumption  that  m 
is  monotone-decreasing  (i.e.,  that  G  is  convex)  means  that  a  deformation  which  occurred 
in  the  distant  past  has  less  influence  on  the  present  stress  than  one  which  occurred  in  the 
more  recent  past,  i.e.  the  memory  fades  with  time. 


Substitution  of  (1.7)  into  balance  of  linear  momentum  (1.2)  yields 

Utt(x,0  =  0uxz(x,t)  -T  f  m(t  -  r)(vxx(x,t)  -  uxx(z.r))dr  +  f(x,t)  (1.13) 

J  —  OO 

or 

Utt(x,t)  =  Euxx(x,t)  ~  f  m(t  -  t)uX£(x,t)(It  +  f(x,t),  (1.14) 

J  —  OC 

where  the  second  form  is  valid  only  if  m  is  integrable.  We  assume  that  the  history  of  u 
prior  to  time  t  =  0  as  well  as  the  values  of  u(x, 0)  and  tit(x,0)  are  prescribed. 

In  much  of  the  work  on  viscoelasticity  it  is  assumed  that  m  is  smooth  on  [0,  oo)  and 
hence  that  m(0)  is  finite.  There  are,  however,  theoretical  and  experimental  indications 
that  for  certain  viscoelastic  materials  it  is  reasonable  to  assume 

m(0)  =  +oo  (1.15) 

(cf.,  e.g.  Doi  and  Edwards  (1978,  1979),  Joseph,  Riccius,  and  Arney  (1986),  Laun  (1978), 
Rouse  (1953),  and  Zimm  (1956)).  Memory  functions  satisfying  (1.15)  will  be  called  sin¬ 
gular.  We  always  assume,  without  explicit  mention,  that  m  is  smooth  on  (0,  oo)  and  that 
(1.8)  holds.  However,  we  shall  always  state  explicitly  whether  or  not  m  is  required  to  be 
smooth  on  [0, oo).  For  singular  memory  functions,  an  important  distinction  arises  as  to 
whether  or  not  the  singularity  is  integrable,  i.e.  whether  or  not  m  is  integrable  at  zero. 
We  note  that  if  m  has  a  nonintegrable  singularity,  e.g.  m(t)  =  e~it~3^2,  then  G(0)  =  +  oo. 

Singular  memory  functions  have  been  employed  in  applications  to  many  important 
engineering  problems  (cf.,  e.g.  Bert  (1973)  and  Walton,  Nachman  and  Schapery  (1978)). 
The  possibility  of  a  singular  memory  function  was  already  considered  by  Boltzmann  (1876) 
who  proposed  that  m(s)  ~  s~x  (except  for  large  s)  and  fitted  such  a  kernel  to  torsion 
measurements  on  a  filament  of  glass. 

Nonlinear  viscoelastic  materials  are  characterized  by  constitutive  equations  which  ex¬ 
press  the  stress  as  a  (nonlinear)  functional  of  the  temporal  history  of  the  strain.  We 
consider  only  materials  which  are  simple  (cf.  Noll  (1958))  in  the  sense  that  the  stress  at  a 
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material  point  depends  only  on  the  history  of  the  strain  at  that  same  material  point  (and 
not  on  the  histories  of  the  strain  at  other  material  points  or  on  spatial  derivatives  of  the 
strain).  We  suppress  the  variable  x  in  our  discussion  of  constitutive  equations. 

General  constitutive  theories  for  nonlinear  viscoelastic  behavior  which  do  not  assume 
a  specific  form  for  the  stress  functional  have  been  given  by  Coleman  and  Noll  (1960), 
Coleman  and  Mizel  (1967,  1968),  Saut  and  Joseph  (1982)  and  Wang  (1965).  We  shall 
restrict  our  attention  to  constitutive  relations  of  singie-integral  type  with  one  memory 
function: 

rt 

(1.16) 


'(()  =  F(c(t))  +  f  m(l  -  r)M(c(t),t(T))dr. 
J  —  oo 


Here  F  and  M  are  given  smooth  functions  with  F(0)  =  0,  M(p,  p)  =  0  for  all  p,  and  m  is 
as  above.  The  function  F  measures  the  stress  due  to  a  strain  history  that  is  constant  in 
time;  its  derivative  F'  is  called  the  equilibrium  stress  modulus.  For  viscoelastic  solids  it  is 
natural  to  assume  that  /’/(e)  >  0,  at  least  for  e  near  zero. 

If  m  is  integrable  we  may  rewrite  (1.16)  in  the  form 


r {t)=f  m(£  -  r)//(e(t),e(7-))dr, 

J  —  oo 


(1.17) 


where 


H(e,p)  =  M(e,p )  +  F{e)(  f  m(s)ds)-1. 

J  c 


(1.18) 


If  we  fix  t  and  hold  the  values  of  e(r)  fixed  for  — oo  <  t  <  t,  but  vary  the  present  value 
e(£),  then  the  integral  in  (1.17)  can  be  regarded  as  a  function  of  e(<).  The  derivative  of 
this  function  is  called  the  instantaneous  elastic  modulus;  it  measures  the  change  in  stress 
due  to  an  instantaneous  change  in  strain.  We  note  that  the  instantaneous  elastic  modulus 
depends  on  the  present  value  of  the  strain  as  well  as  the  entire  past  history  of  the  strain. 
The  equilibrium  stress  modulus  for  (1.17)  is  given  by 


[f/,i  (e,<)  +  H,2  (f,<);  f  m{s)ds , 

Jo 


(1.19) 
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where  H, i  and  H, 2  denote  the  derivative  of  H  with  respect  to  its  first  and  second  argument, 
respectively. 

The  special  case  of  (1.17)  that  occurs  when 


ff(t.p)  =  /(')  +  s(p) 


(1.20) 


will  be  considered  in  Sections  5  and  6.  In  order  to  be  consistent  with  the  literature  on  this 
model  we  introduce  the  notations 


V?(f)  =  /(«)  /  m(s)ds,  0(p)  =  -o(p), 
Jo 


(1.21) 


and  write  the  constitutive  relation  in  the  form 


r(*)  -  -  f  m(t  —  r)ip(e{r))dT. 

J  —  OO 


(1.22) 


We  note  that  for  (1.22),  the  instantaneous  elastic  modulus  depends  only  on  the  present 
value  of  the  strain  and  is  given  by  ¥?'(«);  the  equilibrium  stress  modulus  is  given  by 


roo 

<p'[e)  —  V’,(c)  /  m(s)ds. 

Jo 


(1.23) 


For  the  constitutive  relation  (1.17)  the  equation  of  balance  of  linear  momentum  (1.2) 
becomes 

uu{x,t)  =  /  m(t  -  T)H(ux(x,t),ux(x,T))xdr  +  f(x,t),  xeB,  t  >  0.  (1.24) 

J  —OO 

We  assume  that  the  history  of  u  prior  to  time  t  =  0  as  well  as  the  values  of  u( x,0)  and 
ut(x,0)  are  prescribed,  i.e.  we  take  initial  data  of  the  form 


u(x,r)  =  w(x,r).  x  6  B,  r  <  0,  U-25) 

u(x,0)  =  u0[x),  ut(x,0)  =  Ui(x),  (1-26) 


where  w,  u 0,  and  u  1  given  functions.  Unless  specified  otherwise,  we  assume  that  tr,  uq. 
and  tii  are  smooth.  However,  we  do  not  necessarily  assume  that  uL-i  and  ti  1  are  compatible 


with  w  in  the  sense  that  uo(x)  —  u»(x,0_),  u  t(:r)  =  trt(x,0“).  If  the  interval  B  is  finite  (or 
semi-infinite)  then  boundary  conditions  will  be  imposed  at  the  endpoints  of  B\  typically, 
either  the  displacement  or  the  stress  will  be  prescribed  at  each  endpoint. 

Carrying  out  the  differentiation  in  (1.24)  we  obtain 

V 

utt{x,t)  =  (  I  m(t  —  t)H, i  {ux{x,t),ux{x,T))dr)uxx(x,t) 

J  —  OO 

+  /  m{t  -  t)H,2  {ux(x,t),ux(x,T))uxx(x,T)d.T 
+  /(x,t),  xeB,  t>  0.  (1.27) 

Observe  that  the  coefficient  of  uxx(x,t)  is  precisely  the  instantaneous  elastic  modulus. 

The  paper  is  organized  as  follow's.  Section  2  is  concerned  with  linear  wave  propagation. 
Our  purpose  there  is  to  discuss  the  effects  of  memory  on  a  discontinuity  in  the  data,  with 
particular  emphasis  on  the  effects  of  a  singular  kernel.  The  remainder  of  the  paper  is 
devoted  to  nonlinear  problems. 

Section  3  contains  a  brief  discussion  of  acceleration  waves.  Section  4  is  concerned  with 
the  formation  of  singularities  in  solutions  with  smooth  (but  large)  initial  data.  In  Section 
5  we  discuss  the  global  existence  and  asymptotic  behavior  of  classical  solutions  for  smooth 
and  small  data. 

Comparatively  little  is  known  about  the  existence  of  weak  solutions  (i.e.  solutions 
with  shocks);  this  topic  is  currently  under  active  investigation.  Some  remarks  on  these 
matters  will  be  given  in  Section  6. 

2.  Linear  Wave  Propagation 

In  this  section  we  give  a  brief  synopsis  of  some  results  concerning  linear  wave  propaga¬ 
tion.  There  is  a  very  large  literature  on  this  subject  (cf.,  e.g.  the  monograph  of  Christensen 
(1971)  for  references  to  early  work).  We  refer  to  Section  II. 3  of  RHN  for  a  much  more 
extensive  discussion  of  the  results  quoted  here. 
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We  concentrate  on  the  so-called  Rayleigh  problem,  in  which  a  discontinuity  is  intro- 


««(*,*)  =  0uzx  (x,t) 


viscoelastic  body: 

-  uzz(x,  t))<1t,  x  >  0,  t  >  0, 

(21). 

t  >  0, 

(2-1)2 

>  0,  r  <  0, 

(2.1)3 

0,  x  >  0. 

(2.1), 

Observe  that  there  is  a  discontinuity  in  the  data  at  x  —  t  —  0.  The  solution  of  (2.1) 
provides  a  great  deal  of  information  concerning  the  effects  of  memory. 

To  simplify  the  exposition  we  assume  that  m  is  of  class  C°°  (i.e.  infinitely  differen¬ 
tiable)  on  (0,  oo)  and  that 


(— l)fcm(fc)(f)  >  0  Vt  >  0,  k  =  0,1,...,  (2.2) 

where  m ^  is  the  k- th  derivative  of  m.  We  note  that  most  of  the  results  quoted  below 
remain  valid  under  much  weaker  assumptions.  However,  (2.2)  is  very  reasonable  from  the 
viewpoint  of  applications.  Memory  functions  satisfying  this  condition  are  called  completely 
monotonic.  (This  class  of  memory  functions  is  very  popular  in  rheology.)  We  also  assume 
that 

0  >  0,  m  £  0.  (2.3) 

The  behavior  of  the  solution  of  (2.1)  depends  on  the  behavior  of  m  at  zero.  The  case 
when  m  has  an  integrable  singularity  is  especially  interesting  because  smoothing  of  the 
discontinuity  can  coexist  with  finite  speed  of  propagation. 

If  m  is  smooth  on  j0,oc)  then  u  is  smooth  (in  fact  analytic)  for  0  <  x  <  ct,  and  u  =  0 


for  x  >  ct,  where 


(cf.  Berry  (1958)).  Across  the  line  x  =  ct,  u  sustains  a  jump  discontinuity;  the  amplitude 
A  of  the  jump  is  given  by 

A(t)  =  exp!  — m(0)f/2c3/2j  (2.5) 

(cf.  Chu  (1962)).  In  other  words  the  discontinuity  at  x  =  t  =  0  propagates  with  constant 
speed  and  its  amplitude  decays  exponentially.  (For  a  purely  elastic  material  the  disconti¬ 
nuity  would  propagate  with  constant  speed  and  its  amplitude  would  remain  constant.) 

If  we  formally  put  m(0)  =  +oo  in  (2.5)  then  A(t)  becomes  zero  for  t  >  0,  which 
suggests  that  a  singularity  in  m  has  a  smoothing  effect  on  the  solution.  This  is  indeed  the 
case;  the  precise  degree  of  smoothing  depends  crucially  on  the  strength  of  the  singularity 
in  m.  Explicit  examples  of  singular  memory  functions  which  lead  to  varying  degrees  of 
smoothing  in  the  solution  were  given  by  Renardy  (1982)  and  Hrusa  and  Renardy  (1985). 

If  m  is  integrable  then  u  is  analytic  in  the  region  0  <  x  <  cf  and  u  —  0  for  i  >  ct, 
where  c  is  given  by  (2.4).  Pruss  (1987)  has  shown  that  u  is  continuous  across  the  line 
x  —  ct  (and  hence  continuous  on  the  quarter  plane  x  >  0,  t  >  0)  if  and  only  if  m(0)  =  +oo. 
Desch  and  Grimmer  (1988)  showed  that  u  is  of  class  C°°  on  the  quarter  plane  x  >  0,  t  >  0 
if  and  only  if  m  has  a  singularity  which  is  stronger  than  logarithmic.  They  also  showed 
that  if  m  has  a  singularity  which  is  weaker  than  logarithmic  then  u  fails  to  be  of  class  C1 
across  the  line  i  =  ct.  An  example  of  a  memory  function  with  a  logarithmic  singularity  is 
given  in  Hrusa  and  Renardy  (1985);  it  is  shown  that  the  smoothness  of  u  across  the  line 
x  =  ct  increases  in  a  manner  proportional  to  t. 

If  m  has  a  nonintegrable  singularity  at  zero  then  the  solution  of  (2.1)  is  analytic  in 
the  quarter  plane  x  >  0,  t  >  0.  Since  u  is  analytic,  the  speed  of  propagation  is  infinite. 

To  summarize,  if  m  is  smooth  on  [0,  oc)  then  the  discontinuity  at  x  —  t  =  0  persists  for 
all  t  >  0.  but  its  amplitude  decays  exponentially.  A  singularity  in  m  has  a  smoothing  effect 
on  the  discontinuity;  the  degree  of  smoothing  can  range  from  infinitesimal  (u  is  continuous, 
but  not  C1)  all  the  way  to  analytic  smoothing.  The  degree  of  smoothing  increases  with 
the  strength  of  the  singularity  in  m. 


3.  Nonlinearity  and  Acceleration  Waves 


The  results  on  linear  wave  propagation  indicate  that  memory  has  a  damping  effect.  In 
the  absence  of  memory,  a  nonlinear  elastic  response  can  lead  to  the  formation  of  singular¬ 
ities  in  finite  time  from  arbitrarily  smooth  initial  data;  formation  of  singularities  generally 
occurs  even  if  the  data  are  small.  Equation  (1.24)  includes  a  nonlinear  instantaneous 
elastic  response  in  conjunction  with  damping  due  to  memory.  A  great  deal  of  insight  con¬ 
cerning  the  interaction  of  a  nonlinear  instantaneous  response  with  the  damping  effects  of 
memory  is  provided  by  the  work  of  Coleman  and  Gurtin  (1965)  on  the  growth  and  decay 
of  acceleration  w'aves. 

We  consider  equation  (1.24)  with  B  =  IR  and  /  =  0,  i.e. 

utt(x,t)  =  f  m(t  —  T)H(ttz(x,t),ux(z,r))xdr,  x  £  IR,  t  >  0.  (3.1) 

J  —  OO 

(The  analysis  of  Coleman  and  Gurtin  was  carried  out  for  a  much  more  general  class  of 
constitutive  equations.)  Throughout  this  section  we  assume  that  m  is  smooth  on  [0,oo) 
and  that  m  >  0,  H,i  >  0.  We  consider  a  solution  of  (3.1)  that  is  of  class  C2  on  each  side 
of  a  smooth  curve  t  =  A(x),  across  which  u,  uz,  and  ut  are  continuous,  but  the  second 
derivatives  experience  a  jump.  Such  a  singularity  in  the  solution  is  called  an  acceleration 
wave.  We  suppose  that  an  acceleration  wave  is  propagating  into  an  undeformed  medium 
at  rest;  we  assume  that  u[x,t)  =  0  for  t  <  A(x).  Coleman,  Gurtin  and  Herrera  (1965) 
showed  that  A'(x)  =  c-1,  where 

c=(fftl(0,0)  [°°  m(s)ds)1/2 ,  (3.2) 

Jo 

i.e.  the  acceleration  wave  propagates  with  constant  speed.  (This  conclusion  rests  heavily 
on  the  assumption  that  the  wave  is  moving  into  an  undeformed  medium  at  rest.) 

Let  us  denote  by  .4(t)  the  jump  in  across  the  curve  t  =  A(x).  Coleman  and  Gurtin 
showed  that  A  satisfies  the  ordinary  differential  equation 


where 


1  f°° 

a  =  n  (0,°)/  m{s)ds},  (3.4) 

Jo 

0= '  (s-5) 

The  differential  equation  (3.3)  can  be  solved  explicitly.  The  following  two  consequences  of 
(3.3)  are  of  particular  interest;  we  assume  that  a  ^  0  and  j3  >  0: 

(i)  If  |A(0)l  <  jfj  then 

A(t)  — +  0  as  t  — *  oc. 

(ii)  If  |A(0)|  >  and  aA(0)  >  0  then  A(f)  becomes  infinite  in  finite  time. 

The  above  results  suggest  that  the  damping  effect  of  memory  will  dominate  in  equation 
(3.1)  if  the  data  are  small,  while  nonlinearity  in  the  instantaneous  elastic  response  will 
dominate  when  the  data  are  large.  A  number  of  recent  results  along  these  lines  will  be 
discussed  in  the  next  two  sections. 

We  note  that  for  a  purely  elastic  material  H, 2  =  0  and  hence  ~  0.  It  follows  from 
equation  (3.3)  that  A(f)  becomes  infinite  in  finite  time  if  aA(0)  >  0.  Hence  acceleration 
waves  of  arbitrarily  small  initial  amplitude  can  blow  up  in  finite  time. 

An  explicit  example  of  an  acceleration  wave  in  a  nonlinear  viscoelastic  material  was 
given  by  Pipkin  (1966).  Subsequently,  Greenberg  (1967)  established  the  existence  of  solu¬ 
tions  with  steady  acceleration  waves  for  a  very  general  class  of  viscoelastic  materials. 

The  results  of  Coleman  and  Gurtin  require  that  m  is  smooth  on  i0,oo).  In  the  linear 
case  singular  kernels  lead  to  instantaneous  damping  of  discontinuities.  It  is  therefore 
natural  to  expect  that  singular  kernels  may  also  lead  to  a  much  stronger  type  of  damping 
in  nonlinear  problems.  However,  no  results  of  this  nature  have  been  obtained  so  far. 

4.  Formation  of  Singularities 

If  m  is  smooth  on  [0,  oo)  and  m  >  0,  Hti  >0  then  a  relatively  straightforward  iteration 
procedure  can  be  used  to  establish  the  local  (in  time)  existence  of  a  classical  solution  to 


(1.24)  provided  the  given  data  are  sufficiently  smooth.  The  solution  will  exist  up  to  a 
maximal  time  Tmax  <  oo;  if  Tmaz  <  oo  then  second  derivatives  of  u  will  become  infinite 
in  finite  time. 

The  work  of  Coleman  and  Gurtin  (1965)  on  acceleration  waves  suggests  that  if  the 
data  are  sufficiently  large  then  Tmax  will  be  finite,  i.e.  the  local  solution  will  develop 
singularities  in  finite  time.  This  is  indeed  the  case;  several  results  of  this  nature  have  been 
established  over  the  last,  few  years. 

For  purposes  of  comparison,  we  begin  by  reviewing  a  classical  result  of  Lax  (1964) 
concerning  the  quasilinear  wave  equation.  Consider  the  initial-boundary  value  problem 

utt  =  <p(ux)x,  0  <  i  <  1,  t  >  0,  (4.1) x 

u(0,t)  =  u(l,f)  =  0,  t>0,  (4.1)2 

u(x,0)  =  0,  ut(x,0)  =  uj(x),  (4.1)3 

where  p'  >  0,  p "  >  0,  and  u\  is  smooth  and  vanishes  at  the  endpoints  x  =  0, 1.  (We  take 
tx(x,0)  =  0  only  for  the  sake  of  simplicity;  Lax’s  argument  applies  to  general  initial  data.) 
The  problem  (4.1)  has  a  unique  classical  solution  on  a  maximal  time  interval  [0,  Tmax). 
Lax  shows  that  if  ui(x)  ^  0  then  Tmaz  <  00,  i.e.  the  solution  develops  singularities  in  finite 
time  unless  U)  vanishes  identically  (in  which  case  the  solution  also  vanishes  identically). 
As  t  t  Tmax,  the  second  derivatives  of  u  become  infinite.  If  Ui(x)  *  0  then 

A  =  max  u',(x)  >  0  (4.2) 

since  uj(0)  =  ui(l)  =  0.  An  approximation  for  Tmax  which  is  valid  when  ^max^uifx)!  is 
sufficiently  small  is  given  by 

rmaiS=4^,(0)(A^"(0))-1.  (4.3) 

Under  the  above  assumptions  on  p  and  the  boundary  conditions  (4.1)2,  the  result  for 
general  initial  data  is  as  follows:  The  solution  will  develop  singularities  in  finite  time 
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unless  the  initial  data  for  u  and  v <  both  vanish  identically.  (For  general  data  the  analog  of 
(4.3)  is  more  complicated.)  The  situation  is  similar  if  the  assumption  <p"  >  0  is  replaced 


by  <p’!  <  0. 


Lax’s  argument  makes  crucial  use  of  the  fact  that,  for  0  <  x  <  1  and  0  <  t  <  Tmax, 


the  quantity  tp"(ux(x,  <))  never  vanishes.  Results  establishing  formation  of  singularities 
for  equation  (4.1)  i  when  <p  has  inflection  points  have  been  given  by  MacCamy  and  Mizel 
(1967)  and  Klainerman  and  Majaa  (1981). 

Under  appropriate  sign  assumptions,  the  memory  in  (1.24)  induces  a  dissipative  mech¬ 


anism  which  competes  with  the  destabilizing  effects  of  a  nonlinear  instantaneous  elastic 


response.  The  dissipation  dominates  if  the  solution  is  small,  while  the  elastic  response 


dominates  if  the  solution  is  large.  One  can  still  establish  formation  of  singularities  in  finite 


time  for  solutions  of  (1.24),  but  it  must  be  assumed  that  the  data  are  large  enough  so  that 


the  instantaneous  response  dominates.  The  following  result  is  a  special  case  of  a  theorem 


of  Dafermos  (1986). 


Consider  the  problem 


v,tt{x,t)  =  /  m(t  -  t)H(ux{xJ),ux[x,t))x(It,  0  <  x  <  1,  t  >  0, 
J  —OO 


u(0,f)  =  u(l,t)  =  0, 


u(x,t)  =  0,  0  <  x  <  1,  f  <  0, 


u(x,0)=0,  ut(x,0)  =  uj(x), 


where  Uj  is  smooth  and  vanishes  at  the  endpoints  x  =  0, 1.  We  assume  that  m  is  smooth 


on  [0,oo),  m  >  0,  H, i  >  0,  and  H,u  (0,0)  >  0.  The  problem  (4.4)  has  a  unique  classical 


solution  which  exists  up  to  a  maximal  time  Tmax  <  oc.  Dafermos  proved  that  given  any 


numbers  N,T  >  0  there  are  corresponding  numbers  6,M  >  0  such  that  if 


max  |ui(x)l  <  6,  max  u\  (x)  >  A/,  min  u\  (x)  >  -N 

0<x<l  0<  x<  1  o  <  x  < 1 
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then  Tmaz  <  T.  As  t  *  Tmax  the  second  derivatives  of  u  become  infinite.  Results  of  a 
similar  nature  have  been  obtained  by  Kosinski  (1977),  Slemrod  (1978),  Gripenberg  (1982), 
Hattori  (1982),  Nohel  and  Renardy  (1987)  and  Rammaha  (1987);  for  some  numerical  work 
concerning  the  formation  of  singularities  see  Markowich  and  Renardy  (1984). 

All  of  the  results  concerning  formation  of  singularities  for  equation  (4.4) i  require 
that  m(s)  remains  bounded  as  s  j  0.  It  would  be  very  interesting  to  know  whether  or 
not  singularities  can  form  from  smooth  data  for  nonlinear  problems  with  singular  kernels. 
This  question  is  currently  under  investigation. 


5.  Global  Existence 

This  section  is  concerned  with  the  global  existence  and  asymptotic  behavior  (as  t  — + 
oo)  of  smooth  solutions.  We  assume  throughout  that  m  satisfies 


m(t)  >  0,  m'(f)  <0,  Vt  >  0, 


(5-1), 


m  £  0. 


(5.1): 


Much  of  the  work  on  global  existence  has  been  carried  out  for  the  constitutive  relation 
(1.22),  i.e.  for  the  equation  of  motion 

•t 


Utt(x,t)  =  <p(ux(z,t))x  -  f  m(t  -  t)xI)[ux(x,t))xcLt  +  f{x,t). 

J  —  OO 


(5.2) 


We  assume  that  u  vanishes  identically  prior  to  time  t  =  0,  and  we  impose  nontrivial  data 
for  u  and  ut  at  t  =  0.  Thus  we  consider  initial  value  problems  of  the  form 

utt(z,t)  =  <p{ux{x,t))z-  f  m(t-T)tJ;{ux(x,T))xdT  +  f{x,t),  x  €  B,  t  >  0,  (5.3)  i 

Jo 

u(x,0)  =  u0(x),  itt(x,0)  =  U}(x),  x  e  £,  (5.3)2 

together  with  appropriate  boundary  conditions.  We  assume  that 


<  00,  /  tm(t)dt 
Jo 


<  00 


(5.4) 


m(t)dt)il’'( 0)  >  0. 


v>'(0)  >  0,  0'(O)  >  0,  *s'(0)  -  (  f  m{t)dt)xl’'{0)  >  0.  (5.5)2 

Jo 

Note  that  our  assumptions  permit  the  kernel  m  to  have  an  integrable  singularity  at  zero. 

In  order  to  obtain  global  existence  we  require  uo,  ulf  and  /  to  be  small  and  compatible 

\ 

with  the  boundary  conditions.  If  m  is  smooth  on  [0,  oc)  then  the  results  of  the  preceding 
section  indicate  that  we  must  restrict  the  size  cf  the  data  in  order  to  avoid  the  formation 
of  singularities  in  finite  time.  For  m(0)  =  +  oo,  it  is  not  known  whether  one  can  obtain 
global  existence  for  arbitrarily  large  data. 

We  begin  with  the  case  B  =  [0,1].  To  measure  the  size  of  the  data  we  define 

£/(u0,«i)  =  [  [u0(x)2  +  Uo(x)2  +  Uo(x)2  +  Uo'(x)2  +  u1(z)2-i-Ui(x)2  +  Ui'(x)2]di,  (5.6)i 
Jo 

F(f)  —  f  f  {f2  +  fx  +  ft  +  fit  +  /tt}(x,*)<£x<ft.  (5.6)2 

Jo  Jo 

For  the  Dirichlet  boundary  conditions 


u(0,t)  =  u(l,0)  =0,  t  >  0 
the  compatibility  conditions  will  be  satisfied  if  we  assume 


(5.7) 


uo(0)  =  u0(l)  =  ui(0)  =  ux(l)  =  0,  (5.8)i 

«o(0)  =  u£(l)  =  /(0,0)  =  /(l,0)  =  0.  (5.8)2 

We  note  that  if  m  is  smooth  on  [0,  00)  then  (5. 8)2  can  be  weakened  (cf.  Theorem  IV. 5  of 
RHN).  The  existence  theorem  reads  as  follows:  There  exists  a  constant  p  >  0  such  that  if 
«i>  and  /  satisfy  (5.8)  and 

U{u0,ui)  -  F{f)  <  n,  (5.9) 

then  the  initial  boundary  value  problem  (5.3).  (5.7)  (with  B  -  [0,  l])  has  a  unique  smooth 
solution  u  that  exists  for  all  t  >  0.  Moreover,  as  t  —> >  00,  u,  ur,  ut,  uzx,  uzt,  utt  converge 
to  zero  uniformly  in  x. 


A  very  similar  result  holds  for  the  Neumann  conditions 


u*(0,t)  =  «x(U)  =  0,  t>  0. 


(5.10) 


(We  note  that  (5.10)  is  equivalent  to  o(0,t)  =  <r(l,<)  =  0,  t  >  0;  cf.,  e.g.  Section  3 
of  Dafermos  and  Nohel  (1981).)  The  boundary  conditions  (5.10)  permit  nontrivial  rigid 
motions;  in  order  to  eliminate  the  possibility  of  such  motions  we  normalize  the  data  so 


/  uo  [x)dx  =  /  ui(x)dx  =  0 

Jo  Jo 

f  f(x,t)dx  =  0,  Vt  >  0. 

Jo 


(5.11), 


(5.11), 


There  is  no  loss  of  generality  in  assuming  (5.11)  because  a  problem  with  general  data  can 
always  be  converted  to  an  equivalent  problem  in  which  the  data  satisfy  (5.11)  (cf.,  e.g. 
Section  1  of  Dafermos  and  Nohel  (1981)).  In  place  of  (5.8)  we  now  assume 


“o(°)  =  “o(l)  =  “i(0)  =  “'ll1)  =  °. 


(5.12) 


and  the  existence  theorem  reads:  There  is  a  constant  n  >  0  such  that  if  u0,  uh  /  satisfy 
(5.9),  (5.11),  (5.12)  then  the  initial-boundary  value  problem  (5.3),  (5.10)  has  a  unique 
smooth  solution  u  that  exists  for  all  t  >  0.  Moreover,  as  t  — ►  oo,  u,  u*,  ut,  uxx,  uxt,  utt 
converge  to  zero  uniformly  in  x. 

An  analogous  result  holds  if  (5.10)  is  replaced  by 


«*(0,0  =  «(M)  =  °,  t>0, 


(5.13) 


(or  u(0,t)  =  u*(M)  =  0).  Since  the  boundary  conditions  (5.13)  do  not  permit  nontrivial 
rigid  motions  there  is  no  reason  to  require  (5.11).  A  suitable  “mixed”  version  of  (5.8), 
(5.12)  should  be  assumed. 

For  the  case  when  B  =  IR  we  introduce 

/oo 

(uo(z)2  +  Uol1)2  +  «o'(*)2  +  M*)2  +  ui(J)2  +  “l  w2i  dx  (5.14)  1 
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F*(f)  =  (  (/  f{xj)2dxy/2dt)2  +  /  /  {/i  +  /t2  +  /it  +  /«}(*. (5-14): 

J  0  J —oo  JO  J  —  oo 


to  measure  the  data,  and  the  existence  theorem  reads:  There  is  a  constant  \x"  >  0  such 


that  if 


tfo>o,«i)  +  /”(/)  <  M* 


(5.15) 


then  the  initial  value  problem  (5.3)  (with  B  =  IR)  has  a  unique  smooth  solution  u  that 


exists  for  all  t  >0.  Moreover,  as  <  — »  oo,  u2,  u£,  ttZI,  ul£,  ti££  converge  to  zero  uniformly  in 


2.  (We  note  that  even  if  one  adds  u0(x) 2  to  the  integrand  in  (5.14)  j  the  theorem  provides 
no  information  on  the  behavior  of  the  displacement  u(x,t)  as  t.  —*  oo.) 


The  fieorems  quoted  above  are  a  combination  of  the  results  of  Dafermos  and  Nohel 


(1981),  Hrusa  and  Nohel  (1S85),  Hrusa  and  Renardy  (1986),  and  Hrusa  (1988).  The  first 
global  existence  result  for  (5.3)  was  established  by  MacCamy  (1977)  in  the  special  case 


when  rp  =  <p.  Existence  theorems  for  this  special  case  were  also  obtained  by  Dafermos  and 


Nohel  (1979)  and  Staffans  (1980).  Global  estimates  for  (5.2)  with  m(t)  —  exp(-at)  were 
obtained  by  Greenberg  (1977). 

The  results  described  above  can  be  extended  to  the  constitutive  relation  (1.17)  pro¬ 
vided  that  m  satisfies  (5.1)  and  (5.4).  In  place  of  (5.5)  it  should  be  assumed  that 


17(0,0)  =  0, 


(5.16), 


H,i{ 0,0)  >0,  17,2(0,0)  <  0,  i7tl(0,0)  +  tf  2(0,0)  >  0. 


(5.16): 


For  memory  functions  that  are  smooth  on  [0, oo),  a  more  general  class  of  constitutive 


relations  is  discussed  in  Hrusa  (1983,1985).  The  case  of  an  integrable  singularity  in  m  is 


discussed  in  Hrusa  and  Renardy  (1988). 

Renardy  (1988)  has  established  a  global  existence  theorem  for  the  constitutive  relation 
(1.16)  when  m  has  a  nonintegrable  singularity  at  zero.  He  considered  the  problem 


utt(x,i)  =F(ux(z,t))x  +  f  m(t  -  r)M(ux(x,t),ux(x,r))xdr 

J  —  OO 


+  /(x,0,  0  <  I  <  1,  t>0, 


(5.17), 


Tw »\ f  .r 


mii 


-  ‘O  O  O  *>1 


tt(o,t)  =  *(1,0  =  o,  t  >  o, 


u{x,t)  =  w(x,t),  0  <  x  <  1,  t  <  0, 


(5.17); 


(5.17); 


u(i,0+)  =  u;(x,0  ), ut(x,0  +  )  =  wt(x,0  ),  0  <  x  <  1, 


(5.17). 


where  tv  is  a  given  smooth  function  that  satisfies  equation  (5.17)  i  and  the  boundary 
conditions  (5.17)2  for  t  <  0.  Of  course,  one  can  always  make  the  history  w  satisfy  the 
equation  of  motion  (5.17)  i  for  t  <  0  by  defining  f(x,t)  appropriately  for  t  <  0.  However,  it 
is  assumed  that  /  is  smooth  across  t  =  0  and  this  represents  compatibility  of  the  solution 
for  t  >  0  with  the  given  history.  (We  note  that  when  m  has  a  nonintegrable  singularity,  it 
is  essential  to  require  w[x,0~)  =  u(x,0+).) 

The  assumptions  on  F  and  M  are 


F(0)  =  0,  M(p,p)  =  0  Vp  e  IR, 


(5.18), 


-F'(O)  >  = -M,2[p,p)  >  0  Vj>6K. 


(5.18): 


The  memory  function  m  is  assumed  to  satisfy 


m(t)  >  0,  m'[t)  <0  Vt  >  0, 


f  tm(t)dt 
Jo 


<  oo, 


(5.19), 


(5.19); 


and  an  additional  condition  of  frequency  domain  type.  To  state  this  condition  let  us  put 


K{t)  =  f  m(s)ds  V<  >  0, 
r  oo 

K (iw)  =  /  e~lutK{t)dt  Vw  €  IR. 
Jo 


Finally,  we  assume  there  is  a  constant  C  >  0  such  that 


(5.20), 


(5.20); 


\ReK{iu)\  >  C\ImK{iuj)\  Vw  €  IR. 


(5.21) 


(This  assumption  will  be  satisfied  if,  for  example,  rn(t)  ~  t~a  as  t  — *  0+  for  some  o  with 
1  <  a  <  2.)  Under  the  above  assumptions,  there  is  a  constant  u  >  0  such  that  if  /  is 
smooth  on  [0,  ll  x  IR  and 

/  /  {f2  +  f?  +  f?t}{x,t)dxdt  <  v  (5.22) 

J  —  oo  JO 

\ 

then  the  problem  (5.17)  has  a  unique  solution  that  is  smooth  and  small  (in  an  appropriate 
function  space)  and  exists  for  all  t  >  0.  Moreover,  as  t  — >  oo,  u,  uz,  Ut  converge  to  zero 
uniformly  in  x. 

The  assumption  that  /  is  small  (in  the  sense  of  (5.22))  implies  that  the  given  history  w 
must  be  small.  It  is  not  known  whether  global  solutions  exist  for  large  data.  This  problem 
is  currently  under  investigation.  A  local  (in  time)  existence  theorem  for  large  data  is  also 
given  in  Renardy  (1988).  A  related  local  existence  theorem  is  given  in  Hrusa  and  Renardy 
(1988). 

6.  Remarks  on  Weak  Solutions 

By  a  weak  solution  we  mean  a  displacement  field  u  that  is  allowed  to  have  discontinu¬ 
ities  in  the  strain  uz  and  the  velocity  u*,  and  that  satisfies  the  equation  of  motion  in  the 
sense  of  distributions.  Motions  that  are  well  modelled  by  discontinuous  strain  and  velocity 
fields  are  observed  experimentally  (cf.,  e.g.  Nunziato,  Walsh,  Schuler  and  Barker  (1974) 
and  Walsh  (1985)).  Moreover,  in  view  of  the  results  concerning  formation  of  singularities, 
one  must  allow  for  weak  solutions  (at  least  when  m(0)  is  finite)  if  global  existence  for  large 
data  is  desired. 

There  have  been  many  formal  studies  of  shock  waves  in  viscoelastic  media  (cf.  the 
encyclopedia  article  of  Chen  (1973)  and  the  references  cited  therein).  However,  despite 
the  importance  of  the  subject,  comparatively  little  is  known  about  the  existence  of  weak 
solutions. 

For  a  very  general  class  of  constitutive  relations  (with  smooth  memory  functions), 
Greenberg  (1967)  established  the  existence  of  a  special  class  of  weak  solutions  with  steady 


shocks,  i.e.  solutions  of  the  form  u(i.t)  =  g[x  +  ct),  where  c  is  a  constant  and  g'  has  a 
jump  discontinuity.  It  was  pointed  out  in  Section  11.6  of  RHN  that,  under  appropriate 
assumptions  on  H ,  Greenberg's  proof  can  be  applied  to  the  constitutive  relation  (1.17) 
when  m  has  an  integrable  singularity. 

For  memory  functions  with  integrable  singularities  Londen  (1978)  and  Engler  (1987) 
have  established  global  (in  time)  existence  of  weak  solutions.  Their  result  apply  to  very 
general  classes  of  data  (which  permit  the  initial  strain  and  the  initial  velocity  to  be  discon¬ 
tinuous);  no  size  restrictions  on  the  data  are  needed.  However,  growth  conditions  (for  large 
strains)  are  required  for  the  constitutive  functions;  experimental  data  indicate  that  such 
growth  conditions  may  not  always  be  appropriate.  Londen  considers  the  initial-boundary 
value  problem  (5.3),  (5.7)  (with  B  =  [0,  lj)  in  the  special  case  when  xp  —  <P-  Engler  studies 
the  constitutive  relation  (1.17).  He  also  studies  certain  anti-plane  shearing  motions  of  K- 
BKZ  materials  that  can  be  regarded  as  three-dimensional  analogues  of  (1.17);  a  restriction 
on  the  constitutive  relation  is  needed  to  make  such  motions  possible.  The  arguments  of 
Londen  and  Engler  make  essential  use  of  the  assumption  that  m(0)  =  -foo. 

Recently,  Nohel,  Rogers,  and  Tzavaras  (1988)  established  the  (global  in  time)  existence 
of  weak  solutions  to  the  initial  value  problem  (5.3)  with  B  =  1R  in  the  special  case  when 
xp  =  <p.  They  assume  that  m  is  smooth  on  [0,  oo),  <£>'  >  0,  and  that  <p  has  exactly  one 
inflection  point.  Their  theorem  applies  to  a  very  general  class  of  initial  data  (which  permits 
the  initial  strain  and  the  initial  velocity  to  be  discontinuous).  The  proof  makes  essential 
use  of  the  assumptions  that  m  is  smooth  on  [0,  oo)  and  that  xp  =  <p;  the  extension  to  a 
more  general  class  of  constitutive  equations  does  not  appear  to  be  straightforward. 

The  problem  of  developing  a  more  complete  existence  theory  for  weak  solutions  is 
currently  under  active  study.  As  a  preliminary  step,  several  authors  have  studied  weak 
solutions  of  (scalar)  first-order  integrodifferential  equations  (cf.,  e.g.  Greenberg,  Hsiao, 
and  MacCamy  (1981)  and  Dafermos  (1988)). 
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